Coherence lengths for superconductivity in the two-orbital negative-f/ Hubbard model 
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We study the peculiarities of coherency in the superconductivity of two-orbital system. The 
superconducting phase transition is caused here by the on-site intra-orbital attractions (negative-C/ 
Hubbard model) and inter-orbital pair-transfer interaction. The dependencies of critical and non- 
critical correlation lengths on interaction channels and band fillings are analyzed. 
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I. INTRODUCTION 



Multi-component superconductivity models, devel- 

fed more than fifty years starting from the papers [U- 
include rather varied physics. In connection with the 
presence of interacting order parameters the main prop- 
erties of the multi-band systems are quite different from 
the corresponding characteristics in single-band super- 
conductors. The examination of theses peculiarities has 
been an object of growing interest. 

The various multi-component theoretical scenarios 
have been applied for a number of superconducting ma- 
terials (see Refs. [l-Q and references therein). The re- 
search activity in this direction has been especially stimu- 
lated by the acceptance of the multi-gap superconductiv- 
ity in MgB2^, cuprates^ and iron-arsenic compoundsi^. 
In particular, the derivation of high-quality supercon- 
ducting regions from oxygen ordering, observed recently 
in La2Cu04-f.jpi, supports the multi-band theoretical 
scenario of superconductivity in cupratcs. 

In the present contribution we study the coherency of 
the superconducting ordering of a two-band (two-orbital) 
system with the negative-?/ Hubbard intra-orbital pair- 
ing and inter-orbital pair-transfer interaction. One can 
distinguish here two characteristic length scalesi^ in the 
spatial behaviour of superconducting fluctuations. One 
of these lengths as a function of temperature behaves 
critically diverging at the phase transition point. The 
other one remains finite and its temperature dependence 
is weaker. The formation of these length scales is caused 
by the interband interaction mixing the superconducting 
order parameters of initially non-interacting bands. As a 
result the critical and non-critical coherence lengths as- 
sociate with critical and non-critical fluctuations (see e.g. 
Refs. (islfl^ l) which appear as the certain linear combi- 
nations of the deviations from the equilibrium band su- 
perconducting orders. Consequently, these length scales 
cannot be attributed to different bands involvedi^ (see 
also Refs. [iGl - flsj ). Our results have been obtained us- 
ing the superconducting negative-C/ Hubbard modeliS for 
a two-orbital system on a two-dimensional lattice. 



II. TWO-ORBITAL MODEL OF 
SUPERCONDUCTIVITY 

We start with the Hamiltonian of the two-orbital 
model^S of the following form: 

^ = E E E + - ^^,] 

Q i.j (T 

+ ^EEE^""'«--«»'-«--«-'— (1) 



a, a' I 



where a,^^ [o-iaa) is the electron creation (destruction) 
operator in the orbital a = 1,2 localized at the site i; a 
is the spin index; tf" is the hopping integral; £° is the 
orbital energy; /i is the chemical potential; C/"" < is 
the intra-orbital attraction energy; riiaa = af^^ataa is 
the particle number operator; with a ^ a' is the 

inter-orbital interaction energy. 

The transformation of the Hamiltonian (Eq. [ij into 
the reciprocal space leaves us with the expression 

H = X!EE['^°('^)^^]"ak<T""k<T 
a k IT 

^ E E E E aa(k+q)a«"'k- 



2N 

a, a' k,k' q cr 
^Q(k'-q)-o-'^"'k'-cr 



(2) 



where (k) is the electron band energy with a wave vec- 
tor k, associated with the orbital a and N is the number 
of lattice sites (number of atoms). 

For the description of spatially homogeneous super- 
conductivity one introduces the equilibrium mean-field 
Hamiltonian (the terms which do not contain operators 
have been omitted) 

Hmf = EEE^"('')'^ikff""ka 

Q k a 

+ jEEE(^"«ik.aa-k-. + /»-C-) (3) 
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with 

£„(k) =£„(k) + ^C/""71a-M. (4) 

Here the average number of electrons per site 

^ n„ = iV-i ^ ^ ^ (a+k.«"k.)^,,^^ (5) 

a a Ic (7 

determines the position of the chemical potential fi and 
{■■■)h = Z~^Sp...exp{—H/kBT) denotes averaging pro- 
cedure, where Z is a partition function, ks is the Boltz- 
mann constant, and T is temperature. The homogeneous 
equilibrium superconductivity gaps are defined as 

A„ = E (a"'-Haa'kt>H„^^ . (6) 

a' k 

The Eqs. ^ and ([6]) should be solved self-consistently. 
It is easy to find that 
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l-i^l^tanhMl) 
E^{k) 2kBT 



and 



(a.-u,a„kt>«„,=^^tanh^ 



(7) 



(8) 
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FIG. 1: The critical temperature Tc versus band filling n 
(equal for each band n = ni — n2) for the chosen set of 
interactions: [/" = -1.2t, C/^^ = -1.3t. \U^'^\ = \U'^^\ was 
0.01, 0.04, and 0.07 t for 'l'-'3', respectively. 
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where 



(9) 



The system of gap equations ^ may be now presented 
in the form 



with 



C.(T,A„)^(2^)-X:^tanhfM 



(10) 



(11) 



The critical temperature of the phase transition in the 
two-gap superconductor (Eq. 1) under consideration is 
determined by the following condition 







(12) 



0.14 
0.12 
^0.10 
"'0.O8 
0.06 
0.04 
0.02 
0.0 













A2 


Ai\ \ . 


n=07S\ 









(b) 



.00 0.02 



0.04 
T/t 



0.06 0.08 



FIG. 2: Superconducting gap parameters Ai and A2 ver- 
sus temperature T for various band fillings. U^^ = —1.2t, 
= -1.3t. \U^^\ = |[/^' 1=0.01, 0.04 t, for (a) and (b), 
respectively. 



with ga (Tc) = Ca (Tc O)- 

The calculated [23| dependency of superconducting 
phase transition temperature Tc on band filling n is de- 
picted in Fig. 1. As expected, Van Hove singularity 
present in the middle of the band^SiSi is reflected as a 
peak in Tc for half filled system (n = 1) where the chem- 
ical potential passes this singular point. The depen- 
dence of the small interband interaction U^^ is presented 
by three separate lines \U^^\ = \U^^\ = 0.01, 0.04, and 



0.07 t for T'-'3', respectively. The larger \U \ increases 
the critical temperature Tc- 

The superconducting gaps as temperature functions 
are shown in Fig. 2 for various inter-orbital interactions 
and band fillings. Note that for fairly low temperature 
(T << Tc), Ai and A2 are simultaneously scaled by dif- 
ferent intra-orbital interactions Ua (C/n 7^ U22) and the 
density of states around the chemical potential fi (de- 
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pendent on band filling which was here n = 1.0 or 0.75). 
On the other hand, the influence of U^'^ is visible in the 
vicinity of Tc- The stronger J7^^, the larger Ai. This is 
a consequence of the " inter-orbital proximity effect" . 



III. GINZBURG-LANDAU EQUATIONS 

In the non-homogeneous situation one has to retain the 
non-zero momentum of Cooper pairs in the Hamiltonian 
(Eq. [2]) which can be rewritten as 



H 



a k a 

^ 27V ^ ^ ^ ^ "ik<T"a(-k+q)-o- 
a,a' k.k' q cr 

X aa'(_k'+q)-craQ'k'o- • (13) 



The Bogoliubov theorem determines the upper limit for 
the free energy of the system 

Fb - -fcBrin(Z;,^) + {H- Hl^j)^, , (14) 



Z;„^ = Spexp(-^') 



(15) 



by means of the non-equilibrium mean-field Hamiltonian 
= H 51 £«(k)a+kaa„ka 

a k 

+^ E E E ('5«q«ok.<-k+q)-. + h.c) (16) 



a k,q cr 

with non-equilibrium superconductivity gap order pa- 
rameters (5aq- The free energy upper limit Fb can be 
found as an expansion in powers of (J^q and q. After the 
minimization of Fb and certain transformations one can 
obtain the system of equations for the non-homogeneous 
superconductivity gaps AQ,(r) (the Ginzburg-Landau 
equations) in the following form: 



A„(r) = -^[/""' 



E Po.'U'{T) 



d d 



i.i'=i 



dxi dxv 



A..(r), (17) 



where I refers to the Cartesian axis and d is the dimension 
of lattice. In Eqs. (fTT]) 
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2Ar^e,(k) 



£a(k) 
2fcsT 



(18) 



-E 



2VV Y9|A„ 



tanh 



£;„(k) 2kBT 



(19) 



P""' ^ ^ 47V 4- dqi dqv 1 (k) + £a (k - q) 



tanhf^MV tauhf^i^ 
By supposing further that 

Pali' = PaSlV 

one obtains 



(20) 



q=0 



(21) 



A„(r) = -^C/""' 



5a'(T)-:/„,(r)|A„,(r)r 



A„.(r). 



(22) 



IV. COHERENCE LENGTHS 

The small deviations from the bulk values of supercon- 
ductivity gaps 



-qair) = Aa(r) - A^ 
satisfy linearized equations 



(r) = - ^ 5.' (T) -f (T) ^ ^ 



where 



(23) 



?7«'(i"), 



(24) 



(25) 



with Aq, (T) being the solutions of the system of equations 
(Eq. ini). We introduce the coherence lengths ^ as the 
spatial length scales of r]a (r) : 



771^2 (r) ~ exp - 



Ed 
1=1 



(26) 



From Eqs. ((22)) and ([26|) follows the system 
l + f/"5i(T) + (T)r')r?i(r) 
+ (C/i2.g2(T) + C/i2/32(T)r')r72(r) = 
C/2i5i(T) + C/2i/3i(T)r'^m(r) 
+ ( l + [/22.g2(T) + C/22/32(T)r'V?2(r) = 0. (27) 



A„=0 



the condition for the existence of non-trivial solutions of 
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FIG. 3: The soft and rigid coherence lengths (^r, £,s) versus 
temperature T for several sets of parameters: (a) and (c) 
|C/^^| = 0.01; (b) and (d) |C/^^| = 0.04; (a) and (b) n = 1.0, 
(c) and (d) n = 0.75. ^r, £,s are expressed in terms of lattice 
constant a. 



which, 



l + C/".gi(T) + C/"^, [/i2g2(r) + C/i2^ 



0, 



(28) 



determines the coherence lengths of the two-gap super 
conductor under consideration. On the basis of Eq. 
we have two solutions for ^, 



G (T) ± y/G^ (T) - AK (T) 7(r) 
2K (T) 



(29) 



G(T) 



([/-)' [51 (T)/32(r)+. 92 (m(T)] 
[l + C/"5i(T)] C/22/32(T) 

[l + t/2252(T)] [/"/?! (r), 



A'(r) = [l + C/"5i(T)] [1 + ;72252(T)] 



- ([/i2)'5i(r)52(T), 



7(T) 



11tt22 



u 



12\ 



/3i(T)/32(T). 



(30) 



(31) 



(32) 



The results of our specific two-dimensional lattice cal- 
culation of and versus T/Tc are presented in Fig. 
3. The soft or critical coherence length £,s{T) diverges 
at the phase transition point T = Tc, while the rigid 
or non-critical coherence length ^r(T) remains finite (see 
Fig. 3). Furthermore, for fairly small intcr-orbital in- 
teraction |J7^^| = 0.01 one can notice that shows the 
second maximum (see the peak of in Fig. 3a and the 
rise of 5s with lowering T in Fig. 3c) as the memory about 
the lower transition temperature in the intependent or- 
bital (the one with smaller intra-orbital attraction). Such 
a behaviour is absent for stronger inter-orbital interac- 
tion. Additionally, it is worth to notice that the coher- 
ence length is much larger for n — 0.75 that for ri = 1.0 
(sec Figs. 3a, b and 3c, d respectively). This is again 
the effect of the density of states. Evidently, Van Hove 
singularity^ scales strongly Tc and other superconduct- 
ing parameters including the coherence length. 



V. CONCLUSIONS 

We considered two orbital superconductor with a small 
interaction coupling producing possibility of the Copper 
pair transfer between orbitals. The system have been 
modelled by the negative- C/ Hubbard model on the two 
dimensional lattice. The obtained two values of the co- 
herence lengths reflects the multi-orbital mechanism of 
superconductivity and could have some consequence in 
specific spatial properties of such a complex supeconduct- 
ing statsi^iS^ 

Our results show that band filling strongly effects on 
the coherence lengths. Especially, Van Hove singular- 
ity present in the density of states plays an important 
role in scaling the coherence lengths. The tempera- 
ture dependence of the coherence lengths in our lattice 
model confirms the similar calculations in the continuous 
systen*i^. One solution diverges near the phase transi- 
tion point, while the other one is non-critical. The non- 
monotonic temperature dependence of coherence lengths 
is more pronounced if the inter-orbital coupling is suffi- 
ciently weak. 
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